Spin controlled nanomechanics induced by single-electron tunneling 
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We consider dc-electronic transport through a nanowire suspended between normal- and spin- 
polarized metal leads in the presence of an external magnetic field. We show that magnetomotive 
coupling between the electrical current through the nanowire and vibrations of the wire may result 
in self excitation of mechanical vibrations. The self-excitation mechanism is based on correlations 
between the occupancy of the quantized electronic energy levels inside the nanowire and the velocity 
of the nanowire. We derive conditions for the occurrence of the instability and find stable regimes 
of mechanical oscillations. 



The magnetomotive coupling of electronic and mechan- 
ical degrees of freedom induced by an external magnetic 
field has for many years been a standard tool for imple- 
menting electromechanical functionality on the microme- 
ter and nanometer scales [1, 2], Magnetomotive coupling 
relies on the electron charge in the sense that an electri- 
cal current of charged electrons induces a Lorentz force 
on the current carrying conductor, while the motion of 
the conductor itself induces an electromotive force on the 
charged electrons. However, the electron spin is an ad- 
ditional degree of freedom that may influence electron 
transport and thus potentially affect the clcctromcchan- 
ics of a device. In particular, it is well known that Zee- 
man splitting of electronic energy levels affects the cur- 
rent flow in tunnel structures since the electron tunnel- 
ing rates become spin dependent [3-10]. Furthermore, as 
sample sizes are reduced to the nano scale, mesoscopic 
phenomena such as Coulomb blockade of tunneling [11] 
and quantization of electronic energy levels start [12] to 
significantly affect electron transport, leading to a pro- 
nounced nonlinear current- voltage characteristics (CVC) 
[13-15]. In this Letter we demonstrate that the inter- 
play between spintronics and nanomechanics induced by 
such mesoscopic effects gives rise to a fundamentally new 
set of phenomena. In particular, we show below that 
this interplay may result in self excitation of mechani- 
cal vibrations of a suspended nanowire subjected to an 
external magnetic field and a dc voltage bias. 

To be specific we consider the system sketched in 
Fig. la, which shows a single-wall carbon nanotubc 
(CNT) suspended between a normal- and a magnetic 
metal lead. An external magnetic field H is applied per- 
pendicular to the nanotubc and parallel to the magne- 
tization direction in the magnetic lead. We restrict the 
nanotube mechanical degrees of freedom to the funda- 
mental bending mode and model it as a classical har- 
monic oscillator with frequency u>. The time evolution of 
the amplitude of the fundamental bending mode u(t) is 
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FIG. 1: (a) A carbon nanotube (CNT) subject to an external 
magnetic field H , suspended between normal- (n) and mag- 
netic (m) metal leads biased by voltages ±V5,. (b) Electronic 
energy scheme for the junction: fi is the chemical potential, 
£•(■ 4 are spin-split levels in the CNT, and are electron 
energy bands in the leads (k is the wave vector, a is spin). 

given by the equation 

u(t) + ^u{t)+uj 2 u{t) = 9f(r), (1) 
9f(r) = -m- lr Ti(d u H e p{ T )) , 



where Q is the quality factor characterizing the dissipa- 
tive processes induced by the coupling to heat baths in 
the leads, m is the mass of the oscillator, p is the den- 
sity matrix of the electronic subsystem described by the 
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Hamiltonian ff e , defined as 
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Here at and c, arc creation operators for electrons in 
the normal- and spin-polarized metal leads, respectively, 
d\, is the creation operator for electrons in the nanotube, 
while the index a =t, I labels the electronic spin. 

The Hamiltonian Hi in Eq. (2) describes the two leads 
as reservoirs of non-interacting electrons, while the ex- 
pression for Hd describes the electronic properties of the 
nanotube considered as a quantum dot, with its second 
term representing the Coulomb interaction between two 
electrons with different spin projections. The spacing be- 
tween electronic levels inside the nanotube is assumed to 
be large enough with respect to the symmetrically ap- 
plied bias voltage 2V& for only a single pair of Zeeman- 
split electron energy levels, £-(-4 = e Q ± pbH/2, to be 
relevant for electronic transport (if if = the energy 
level eo is spin-degenerate; /is is the Bohr magneton). 

The Hamiltonian H t in Eq. (2) describes electron tun- 
neling between the nanotube and the leads in terms of the 
tunneling amplitudes t n and t m between the nanotube 
and the normal and magnetic lead, respectively. The 
phases of the tunneling amplitudes depend on the nan- 
otube deflection [1G], the parameter pi, = aeHL/c gives 
the nanotube momentum change induced by the Lorentz 
force when an electron tunnels from the nanotube to a 
lead, and a ~ 1 is a numerical factor determined by the 
spacial profile of the fundamental mode. 

The density matrix p obeys the Liouvillc-von Neu- 
mann equation. Introducing the interaction represen- 
tation, p{r) = e -iH T/h~( T yH 0T /h ; H Q = Hi + H d , 
and exploiting the reduced density matrix Ansatz [17] 
p(r) = /Od(r) • pi one gets the following equations for pd 
and the force term, S, in Eq. (1): 
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where Tr; denotes a trace over all electronic states in the 
leads, while pi = e -(Ht-^N)/k B T /Tr^e-^-^/^}, ^ 
is the chemical potential, and N is the electron number 
operator in the leads. 
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The two equations (3) describe the coupled dynamics 
of the electronic and mechanical subsystems. To demon- 
strate the interplay between spintronics and nanome- 
chanics we analyze them for the case when the bias volt- 
age is applied in such a way that the normal lead serves 
as a source (S) electrode while the magnetic lead acts as 
a drain (D) (see Fig. lb). 

Two additional conditions are assumed to be satisfied: 
First, that U 3> eVt ~ e a = s a — p 3> ksT, which means 
that the charge transfer from the drain electrode to the 
nanotube is exponentially suppressed, while the strong 
electron-electron interaction inside the nanotube imposes 
the Coulomb blockade constraint Tr(dl^d^d^d^pd) = 0. 
Secondly, that w <C ^s{D), hsT/h. Here ^s(D) = 
2Tih~ Y t 2 n ^ m ^v^ m y and the electronic density of states 
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are assumed to be independent of energy for both 
leads and also independent of spin in the source lead, so 
that v\ = and hence Tg = Ig . 

In the case of a low mechanical oscillator frequency one 
can use a quasi-static approximation to solve Eqs. (1) - 
(3). Neglecting the time dependence of pd and using the 
approximation u(t) — it(r') ~ u(r)(r — r') one finds that 
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To arrive at Eq. (5) we assumed the magnetic lead to be 
fully spin polarized, so that vj n =0 [13-15]. To assume 
partial polarization would be more realistic, but not im- 
portant for our conclusions, as will be discussed below. 

The physical interpretation of this result is obvious: 
the force acting on the nanotube is the Lorentz force in- 
duced by the quasi-stationary current J(V), while the 
effective voltage V is comprised of the bias voltage Vb 
and the electromotive force aHLii that is induced by the 
motion of the nanowire in the external magnetic field. 
Linearizing Eq. (1) with 3 given by Eq. (4) one finds that 
the electromotive coupling induces an effective damping 
or pumping of mechanical vibrations in accordance with 
the sign of the differential conductance J = dJ/dV. 
In the case of pumping, which corresponds to J < 0, 
one can expect an electromechanical instability resulting 
in self-excited oscillations of the nanowire if Q is large 
enough and the instability parameter j3 is positive, i.e. if 
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> 0. 
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It is obvious from Eq. (5) that the differential conduc- 
tance J (V) is negative for voltages such that \V— £f/e| < 
ksT/e. Even if the magnetic lead is less than 100% po- 
larized, as assumed in (5), J (V) would still be negative 
in the same voltage interval [13-15]. Therefore, the nec- 
essary condition J (V) < for an instability to occur 
does not require 100% polarization. Nevertheless, for the 



sake of clarity, we assume this to be the case in what 
follows. 

The mechanism of the instability is most transparent 
when eVb = e-f, as in Fig. lb, and psH 3> fc^T. Under 
such conditions electrons that jump from the S-lead to 
the spin-up level in the nanotube can only jump back, 
since they can not reach any spin-up state in the D-lead. 
Such processes lead to random changes of the nanotube 
momentum by ±pl , which results in a diffusive motion in 
its phase space. However, an electron that tunnels to the 
spin-down level under the same circumstances can only 
proceed by tunneling to the D-lead. This is because there 
are no empty states in the S-lead at accessible energies. 
Tunneling through the spin-down channel changes the 
net nanotube momentum by pi and hence its mechan- 
ical energy by plu(t). Therefore, correlations between 
the probability for tunneling through the spin-down level 
and the direction of the nanotube motion can give rise to 
cither damping or pumping of its mechanical vibrations. 
Precisely such correlations appear due to the Coulomb 
blockade of tunneling to the spin-down level when the 
spin-up level is occupied. Indeed, the population of the 
spin-up level increases when the nanotube moves in such 
a direction that the electromotive force shifts the "effec- 
tive" chemical potential, ^ ff = p — aHLii, up. A shift in 
this direction means that electron tunneling to the spin- 
down level is decreased due to the Coulomb blockade. 
When the nanotube moves in the opposite direction the 
situation is reversed and the probability for electrons to 
tunnel to the spin-down level is enhanced. This means 
that, on the average, momentum transfer to the nanotube 
is more probable if the nanowire moves in this particu- 
lar direction. In short, the spin- up level in the nanotube 
serves as an effective gate, which — depending on the 
direction in which the nanotube moves — increases or 
decreases the probability for electron tunneling from the 
source- to the drain electrode via the spin-down level; the 
only available channel in our case. The described correla- 
tions between the tunneling processes and the nanotube 
velocity lead to a "pumping" of energy into the mechan- 
ical subsystem. If the rate of pumping is larger than 
the rate of dissipation (due to damping) , the result is a 
mechanical instability. 

To analyze how the instability described above evolves 
we impose the Ansatz u(t) = uq + A(t)cos(u)t) . Here uq 
is the stationary shift of the nanowire position induced 
by the average Lorentz force and A(t) is the amplitude of 
the harmonic oscillations, which by assumption is slowly 
varying, so that A <C wA. Inserting this Ansatz into 
Eq. (1), using Eqs. (4) and (5) and averaging over the 
fast oscillations one gets an equation for the amplitude 
of the form 
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FIG. 2: Behavior of the function defined in Eq. (7) 

for different magnetic fields leading to (a) soft excitation of 
nanowire oscillations at AV = Vb — (eo — p)/e = 0.34 mV 
and (b) hard excitation at AV = 0.41 mV. Other parameter 
values are T = 0.2 K, Q = 3 ■ 10 4 , F l D /Ts = 0.4, u/T s = 0.2 



The first term in this expression for $(A) describes the 
damping of mechanical oscillations due to a coupling to 
the thermal environment, while the second term describes 
the pumping generated by the electronic current. If the 
condition (6) is satisfied, pumping dominates over damp- 
ing at small vibration amplitudes A and hence the am- 
plitude increases with time according to Eq. (7). How- 
ever, as can be inferred from Eq. (5), at large amplitudes 
A 3> ksT/eHLui the pumping term in Eq. (7) saturates 
at ~ (HL/moj)Jo, where Jo ~ e ^r> is the characteristic 
current through the system. Accordingly, the damping 
and pumping terms cancel each other at a finite ampli- 
tude A st ~ (QHL/mu; 2 )J . 

The exact value of the amplitude of stationary oscilla- 
tions can be obtained by solving the equation & (A st ) = 0. 
A numerical analysis of the function &(A) shows that, 
depending on the bias voltage, the onset of stationary 
nanowire oscillations in an increasing magnetic field can 
be either soft (A st increases continuously from zero) or 
hard (A st jumps to a finite value), see Fig. 2. In order 
to investigate this situation analytically we assume the 
instability parameter to be small, \/3(Vb,H)\ <C 1, and 
expand Q(A) in a Taylor series. Keeping terms up to the 
third order in A one finds 
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a) Soft excitation of nanowire oscillations, Fig. 2a. It is 
readily seen from Eq. (8) that the onset of nanowire 
oscillations is soft in the interval V* < Vb < V% , 
where J (Vb) > (here V^ 2 are the points of inflec- 
tion of the negative differential resistance curve J (V)); 
for f3(Vb,H) > the non-moving state of the wire 
is unstable and stationary oscillations with amplitude 
Ag t oc \[Wl 6 " appear spontaneously, smoothly increas- 
ing in amplitude with increasing values of P(Vb, H) (that 
is with increasing H). 
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FIG. 3: (a) Stability diagram for nanowire oscillations in the 
(H, AV) plane, where AV = Vb — (eo — ^) /e. In the white and 
dark grey regions there is only one stable stationary state cor- 
responding to a non-moving and a vibrating nanowire, respec- 
tively. The light gray region is where both these states are sta- 
ble, (b) Average current in units of eFs, in (b-1) as a function 
of bias voltage AV for H — 10 T [horizontal line in (a)] and in 
(b-2) as a function of magnetic field H for AV = 0.4 mV [ver- 
tical line in (a)]. The full (dashed) lines show the result for 
an ascending (descending) magnetic field. The gray dotted 
curves correspond to a static nanowire. Other parameters: 



T = 0.2 K, 



3 ■ io 4 , i%/r s = o.4, oj/Vs = o.i. 



b) Hard excitation of nanowire oscillations, Fig. 2b. 
With bias voltages V& for which J (Vb) < (i.e., Vb < 
V* ,Vb > V 2 * ) , the instability develops in a qualitatively 
different way. There are now two bifurcation points: 
H^\Vb), at which the curve ®(A,H) "touches" the A- 



axis from above (see Fig. 2b) and h[ 2) {V b ) > H { c L '(V b ) 
at which f3(H^ 2) (Vb)) = 0. As long as H < H { c 1] \V b ) the 
non-moving state of the nanowire is stable, as shown in 
Fig. 3. However, in the interval H^\v b ) < H < H^ ] {V b ) 
the system is bistable since here both the non-moving 
state (A = 0) and nanowire vibrations of a finite ampli- 
tude A st ~ (QHL/moj 2 )Jo (see above) are stable. When 
H > m'\V b ) the non-moving state is unstable and vi- 
brations with a finite amplitude is the only stable sta- 
tionary state of the nanowire. 

The existence of a region of magnetic fields and bias 
voltages for which the system is bistable with respect to 
the vibration amplitude results in a hysteretic behavior of 
the averaged current J = J J{V b + aHLujAsm(j))d<fi/2ir 
under a change of magnetic field or bias voltage. Figure 3 



shows that the width of the hysteresis loop is ~ 0.1 mV 
with respect to the bias voltage and ~ 1 T with respect 
to the magnetic field, while its height is ~ J Q . 

In conclusion, we have considered the mechanical prop- 
erties of a nanowire suspended between two metal leads, 
one nonmagnetic and one magnetic. We have shown 
that interplay between Coulomb blockade of tunneling 
and spin-dependent single-electron tunneling gives rise 
to a fundamentally new response to a magnetomotive 
coupling of the electrical current through the nanowire 
and mechanical oscillations of the nanowire. In partic- 
ular, we have demonstrated that in the presence of an 
external static magnetic field, mechanical vibrations of a 
dc voltage-biased nanowire can be self-excited. In con- 
trast to the resonant excitation of nanowire oscillations 
that may be induced by an electrical or magnetic ac sig- 
nal, the amplitude of these self-excited stationary vibra- 
tions is not limited by any resonant condition and, as 
a result, they can be large. In a realistic experimental 
situation with a CNT resonator of length L ~ 1 ^m, vi- 
bration frequency u ~ 2tt x 200 MHz, and quality factor 
Q <~ 10 4 carrying a characteristic current Jo ~ 1 nA, the 
vibration amplitude is A st ~ 10 nm at T ~ 0.1 K in a 
magnetic field of ~ 10 T. In principle, oscillations with 
such an amplitude could be directly monitored by clever 
imaging techniques [18]. We have also demonstrated that 
the onset and disappearance of these mechanical vibra- 
tions are manifest in a pronounced hysteretic behavior 
of the averaged electrical current through the structure. 
This hysteresis also facilitates detection of the self-excited 
nanowire vibrations in an experiment. In general, spin 
accumulation at the nanowire-magnetic lead contact may 
also result in an instability or hysteretic behavior of the 
current-voltage characteristic due to the spin transfer 
torque effect [19]. However, this requires [20] current 
densities ~ 10 /iA/(nm) 2 , which are orders of magnitude 
higher than those considered here. 
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